Modified Li chains as atomic switches 
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Abstract 

We present electronic structure and transport calculations for hydrogen and lithium chains, using 
density functional theory and scattering theory on the Green's function level, to systematically 
study impurity effects on the transmission coefficient. To this end we address various impurity 
configurations. Tight-binding results allow us to interpret our the findings. We analyze under 
which circumstances impurities lead to level splitting and/or can be used to switch between metallic 
and insulating states. We also address the effects of strongly electronegative impurities. 
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I. INTRODUCTION 



Electronic devices have been reduced more and more in size over the last decades. Fur- 
thermore it is now possible to place atoms or molecules accurately between macroscopic 



electrodes, hence experimenta 



studies of the electronic transport for single atoms 



molecules 13*1 , and nanowires 



0,3, 



4j have become available. One-dimensional structures are 
of particular interest due to their restricted transport channels, making them prototypical 
model systems. For example, mono-atomic chains have been realized by molecular beam 
epitaxy Using the tip of a scanning tunneling microscope, it has been possible to place 
a row of eight atoms on a NiAl substrate j^] . 



From the theoretical point of view, transport through distorted one- dimensional systems 
has been addressed already in the early 1990ies by bosonization techniques js-8]. Transport 
properties of lattice models are currently investigated on several levels. Comparison of an 
exact treatment by the density matrix renormalization group [9j] with approaches using den- 
sity functional theory 10] shows that the latter approach often is sufficient to calculate the 
linear conductance, at least qual itatively. The generalization of density functional theory 
to time- dependent potentials [TjJ allows to study the propagation in time of the electronic 



states 



12J. To model the experimental setup more realistically, several methods have been 



developed for describing the transmission through nano-contacts. Tight-binding formula- 

nn 

tions have been applied to metallic nano-contacts |13l. 1141] . adding orbital information from 
chemical analysis; however, these may fail in the contact regime. Nowadays, most appo- 
raches rely on a combination of density functional theory and a scattering approach on the 



non-equilibrium Green's function level, based on the Landauer-Biittiker scheme 



15|- 



In this article, the electronic structure and the transmission coefficient are determined for 
H and Li chains with defects. In Sec. II we present details of the calculational method and 
the structural setup, and discuss in Sec. Ill some results for H chains and their dependence 
on the calculational parameters. In Sec. IV we turn to impurities in Li chains, focusing on 
the interrelation between the energy level spectrum of the scattering region and the transport 
properties. 
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II. COMPUTATIONAL DETAILS AND STRUCTURAL SETUP 



We calculate the transmission through nanowires using a combination of scatte ring theory 



and density functional theory 



[l6, Q, based on the SMEAGOL and SIESTA 



1* 



codes. 



A single zeta basis set and the local density approximation for the exchange-correlation 
potential are used. In order to model the experimental situation (without gate voltage), 
metallic leads are connected to a central scattering region. In our case, the leads are H 
chains, see Fig. [TJ In the Landauer-Biittiker formalism, the self-energies of the left (L) 
and right (R) leads are calculated first. The screening within the metallic leads ensures that 
effects of the contact region decay within a few nanometers. Since the leads are connected to 
the central scattering region (molecule, nano-contact, or interface), an effective description of 
the central region (C) emerges which includes the properties of the leads. In linear response, 
the transmission coefficient is given by the retarded Green's function Gc of the central region 
and the lead self energies T> L / R . With T L / R = i[E L / R (E) — Yj l , r (E)\ we have [3] 

T(E,V = 0)=Tr[T L GlT R Gc}. (1) 
The conductance is given by G = ^-T(Ep), where the factor 2 accounts for the spin 



degeneracy [15 ]. 



The impurity models studied in Sec. IV are displayed in Fig. [TJ The parent structure is 
a finite (metallic) Li chain with three atoms, which is coupled to H-chain leads. By varying 
the distance between the last H and the first Li atom, we change the coupling from strong 
(a = 2.8 A) to weak (a = 4.0 A). The distance between the Li atoms within the chain is 
fixed to b = 3 A. The impurities are Li atoms adjacent to the parent chain, thus breaking 
the rotational and/or inversion symmetry of the system. In configuration (a) an additional 
Li atom is placed directly below a chain atom, compare {20)]. The distance d between the ad- 
atom and the chain is varied. A second atom on the other side of the chain, configuration 
(b), restores the symmetry. An additional atom can also be placed below a Li-Li bond, 
configuration (c), which also breaks the symmetry in the transport direction. In addition, 
we study configuration (d) in which we place the ad-atom below the central bond, to obtain 
a system with inversion symmetry. 

As the transmission coefficient shows resonance peaks at the energy levels of the scattering 
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FIG. 1. Ad-atom configurations under consideration. The distance d between the ad-atom(s) and 
the chain is varied. The i-parameters refer to the tight-binding description of the central scattering 
region; see main text. 



region, an analysis of this spectrum allows an explanation of the behavior of T(E). For this 
purpose, the energy levels can be discussed on the basis of a tight-binding approximation with 
nearest-neighbor hopping. With the hopping parameters indicated in Fig. [I], the Hamiltonian 
for a three-atom Li chain is given by 

H = -t par (c+c 2 + c+c 3 + h.c.) , (2) 

where the sites are numbered 1, 2, and 3. The hopping energy along the chain, £ par , is 
of the order of 2 eV, as the bandwidth of a Li chain with atomic spacing b = 3 A is 
~ 4.5 eV. The three-atom chain has the energy levels £^3 = ±-v/2t par , and €2 = 0. The 
coupling to the leads can be taken into account by an additional hopping to the left and 
to the right of the scattering region. The hopping to the ad-atoms is denoted by t peT , thus 
Hper = — £pcr(c^~C4 + h.c.) for configuration (a), and H per = — t per (c^C4 + c^c^ + h.c.) for (b). 
In case of configurations (c) and (d), the hopping terms are given by — £ s (cf c^ + c^ C4 + h.c), 
and — t s (c^c 5 + C3 c 5 + h.c), respectively. 

4 



III. HYDROGEN CHAINS 



We start our discussion with a review of results for H chains, which is the simplest case 
of one- dimensional scatterers, see the structural setup in Fig. [2J He and coworkers [21] 
have used a system of infinite H chains as leads and a finite six-atom H chain for their 
transport calculations, focusing on the influence of the exchange-correlation potential and 
the according self-interaction errors. In the following we use the local density approximation, 
since it reproduces well the qualitative behavior of the transmission coefficient 10( . Clearly, 
T(E) depends on the energy levels of the central H chain, and on the coupling strengths to 
the leads. The number of atoms within the chain determines the number of energy levels, 
and hence the number of resonances, see Fig. El A five-atom chain is found to be metallic, 
whereas a six-atom chain turns out to be insulating. This odd-even behavior is reflected 
by a maximum or minimum in T(Ep) in Fig. |2j Moreover, the transmission peaks become 
broader with increasing coupling to the leads, and also shift in energy. 

Considering in addition the bond length within the chain, b, as a parameter, we find, as 
expected, a decrease of the band width with increasing b. For b = 1.0 A the bands are close 
to parabolic, and the unoccupied levels are more sensitive when decreasing the coupling than 
the occupied levels. With increasing b, the dispersion becomes more and more symmetric 



22] . Keeping the lead-chain 



and cosine shaped, consistent with the tight-binding description 
distance fixed at, say, a — 1.2 A but increasing the bond length b beyond this value, one 
arrives at a situation where the first and the last atom of the central region are effectively 
bound to the leads, i.e., the scattering region is a diluted chain, however, with two atoms 
less I22I, |23]. 



IV. LITHIUM CHAINS 

When a finite Li chain is coupled to H chain leads, charge is transfered between the leads 
and the scattering region. For an even number of Li atoms this charge transfer results in 
a level shift towards E F j2l|. Extending the studies of 21], we discuss in the following the 



influence of perturbations. For configuration (a), with an additional Li atom adjacent to the 



Li chain, the tight-binding energy levels of the scatterer are given by £i )4 = ±^2tp ar + t^ er , 
and e 2 3 = 0. The coupling to the leads lifts the degeneracy of the zero-energy level. A 





E - E F (eV) E - E F (eV) 



FIG. 2. Transmission coefficients of five-atom (left) and six-atom (right) H chains coupled to H 
leads, upon increasing the distance a between the scattering region and the leads. The bond lengths 
in the scattering region and in the leads are fixed at 1.0 A. 

comparison of the density of states (DOS) obtained by SIESTA with T(E) is given in Fig. 
|3|a) for a projection of the DOS onto the Li atoms. The distance of the ad-atom to the 
chain is d — 3.5 A. The four energy levels of the scatterer are located at —1.5 eV, eV 
= Ep, 0.3 eV, and 1.7 eV. The smaller peak at —1.0 eV as well as the shoulder at 1.6 eV 
arise from hybridization with atoms in the leads. When the distance between the ad-atom 
and the Li chain is increased, the splitting of the states near Ep becomes more pronounced, 
see Fig. [3](b). 

Analyzing the spatial distribution of the charge for each energy level yields strong similar- 
ities between the tight-binding model and the DFT calculation. The lowest level corresponds 
to a homogeneous charge distribution over the scatterer. The second level (highest occupied 
molecular orbital, HOMO) resides on Li atoms 1, 3, and 4, whereas the third level (lowest 
unoccupied molecular orbital, LUMO) is dominated by the ad-atom, see Fig. Etc). In the 
range d = 3.5 . . . 4.0 A a redistribution of weight from the HOMO to the LUMO is apparent. 
For increasing distance d, the charge becomes more localized on the ad-atom, resulting in 
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FIG. 3. (a) DOS (divided by 4, for easy comparison, in units of 1/eV) and transmission coefficient 
of configuration (a) for d = 3.5 A. (b) Transmission coefficient for various distances between chain 
and impurity, (c) Charge density iso-surfaces of the HOMO and the LUMO. (d) Zoom of (b) 
around Ef- 



a stronger transmission through the LUMO. The minimum of the total energy appears for 
d = 2.75 A; then the energy increases almost linearly when increasing d (up to d — 4.0 A). 

For the symmetrical configuration (b), see Fig. HJ we find five energy levels in the tight- 
binding spectrum: bonding and antibonding levels £i j5 = ±^/2tp ar + 2t^ cr as well as £2,3,4 = 
0. A next-nearest-neighbor hopping t nn shifts one of the zero-energy levels, say, £4, to a finite 
value, given by e 4 ~ 2t nn tp er /t par (which holds in the limit t n n,^ P er "C £ par )- The coupling to 
the leads also lifts the threefold degeneracy of £2,3,4- As a result, a spectrum of five levels 
is obtained, with two nearly degenerate levels near the Fermi energy, see Fig. HI Level 4, 
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E - E F (eV) E - E F (eV) 

FIG. 4. Transmission coefficient for configuration (b) with strong (a = 2.8 A, left) and weak 
(a = 4.0 A, right) coupling to the leads. 

corresponding to £4, has a rather homogeneous charge distribution over the scattering region 
for a small distance, d = 2.5 A. With increasing distance, the contribution of atom 2 to the 
charge density decreases, and hence the transmission through this level. 

The position of the fourth peak in T(E) depends strongly on the distance between the 
ad-atoms and the chain. For d = 3.5 A, for example, it is located at about 0.5 eV, whereas 
it is found at 1.0 eV for d = 3.0 A (equilibrium distance). The two peaks near Ep cannot be 
resolved on the left hand side of Fig. H] due to the coupling to the leads. For weak coupling 
a three-peak structure around Ep is resolved, again showing the redistribution of weight 
below and above Ep, compare Fig. El^d) with the right hand side of Fig. |H 

Configuration (c) corresponds to a scatterer which is asymmetric in both in-plane di- 
rections. The tight-binding model yields the spectrum e 14 ps ±y/2t pSLI and £2,3 ~ ±£ s /\/2 
provided t s is small t s <C £ par , i.e., when the distance between ad-atom and chain is large. 
The first guess thus is that the system should be insulating, in contrast to the parent chain 
and to the configuration where the ad-atom is attached to a single chain atom. According to 
Fig. |5l however, T(E) does not show a minimum at Ep: in fact, due to the charge transfer 
with the leads the HOMO is fixed near Ep, as demonstrated in [21] (also within DFT). In 
addition, the lowest level shifts slightly in energy when varying d. 

Note that for d = 2.5 A the distance between ad- atom and chain atoms is only about 
d = 2.9 A hence smaller than the intra-chain nearest-neighbor distance of d = 3.0 A. 
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E - E F (eV) E - E F (eV) 

FIG. 5. Transmission coefficient for configuration (c) with strong (a = 2.8 A, left) and weak 
(a = 4.0 A, right) coupling to the leads. 

Nevertheless, T(E) for d = 2.5 A essentially is the same as the transmission for larger 
distances. For d above 3.0 A the HOMO and its transmission peak remain unchanged; 
however, the LUMO shifts to lower energy with increasing d (i.e., for decreasing t s ). 

Symmetry in the transport direction can be obtained by placing the ad-atom below 
the central bond of a four-atom Li chain, see configuration (d). In this case the parent 
chain is insulating as £1,2,3,4 = (±1/2 ± \/5/2)tp ar . Including the impurity yields = 
— ipar/2 ± A/5tp ar + 8tg/2, e 2 ,5 = t pax /2 ± v / 5i P ar/2, and e 3 = 0. Thus, the ad-atom again 
switches the electronic state of the chain, but now from insulating to metallic. If we couple 
the central chain to the leads, £3 is shifted to higher energy, due to charge transfer, see Fig. 
|6j For d > 3.0 A the energy of the LUMO and its transmission peak remain essentially 
unaltered, but the distance to the HOMO decreases as t s decreases, similar to the behavior 
of configuration (c). 

We next study the influence of an impurity different from the atoms in the central region, 
comparing ad-atom and intra-chain configurations. As a prototypical example with high 
electronegativity we choose a fluor atom, which we place either adjacent to or incorporated 
into the Li chain. In the former case we expect, for an F atom attached to a single chain 
atom, that one Li electron will be trapped in the Li-F cluster, no longer contributing to 
transport. According to the left hand side of Fig. [7J the transmission spectrum of a seven- 
atom Li chain with F ad-atom in the middle (d = 2 A away from the chain) consists of 
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E - E F (eV) E - E F (eV) 

FIG. 6. Transmission coefficient for configuration (d) with strong (a = 2.8 A, left) and weak 
(a = 4.0 A, right) coupling to the leads. 




E - E F (eV) E - E F (eV) 

FIG. 7. Left: Transmission coefficient for a seven-atom Li chain with an F ad-atom attached 
directly to the central Li atom, and for a six-atom Li chain with weak link at the center. Right: 
Transmission coefficient for a Li chain with intra-chain F impurity. 

six peaks, confirming the above picture of one trapped electron. We note that the chain 
remains metallic despite the impurity. This can be understood by studying an equivalent 
system with a six-atom Li chain and a weak link at the center. To model the weak link, we 
increase the distance between the two central Li atoms to 5 A. Effectively, we obtain the 
transmission spectrum of two three-atom Li chains, see the dashed curve on the left hand 
side of Fig. [TJ i.e., metallic behavior and a high transmission at Ep. 
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The situation is different when the electronegative atom is placed within the chain, as 
shown on the right hand side of Fig. [7] for a Li-Li- F-Li configuration with a Li-F distance 
of 1.5 A. We find that the F atom forms a site with large attractive potential, since one 
transmission peak appears at very low energy (at about —8.2 eV). We note that the resonance 
level at Ep does not attain a transmission of 1. While an F ad-atom forms a neutral cluster 
with the adjacent Li atom, an intra-chain F atom couples to both Li neighbors and acts as 
strong attractive potential. 

V. CONCLUSION 

We have investigated the influence of impurities on the transmission of metallic and 
insulating monovalent Li chains. We have demonstrated that an impurity connected to 
a single atom of the chain leads to a splitting of the zero-energy peak, causing a strong 
dependence of the conductance on structural details of the scatterer. An impurity connected 
to two atoms of the chain with equal distance, on the other hand, switches the system from 
insulating to metallic behavior, or vice versa. The transmission coefficient, however, in 
both cases shows a minimum at Ep as a consequence of the charge transfer with the leads. 
Charge transfer with the leads prohibits a prediction of details of the transmission spectrum 
by tight-binding models of the scatterer. Electronegative ad-atoms localize electrons in the 
chain but do not decouple it, whereas electronegative impurities within the chain result in 
low-lying transmitting levels and, thus, shift the spectrum to higher energy. Our results 
demonstrate that switches based on mono-atomic chains show a high variability of possible 
modifications that can be used to tailor the switching behavior. 

ACKNOWLEDGMENTS 

We gratefully acknowledge discussions with S. Schenk and P. Schwab, and thank the 
Deutsche Forschungsgemeinschaft (TRR 80) for financial support. 



[1] E. Scheer, N. Agrait, J. C. Cuevas, A. L. Yeyati, B. Ludoph, A. Martm-Rodero, G. R. 
Bollinger, J. M. van Ruitenbeek, and C. Urbina, Nature 394, 154 (1998). 

11 



[2] K. Hansen, S. K. Nielsen, M. Brandbyge, E. Laegsgaard, I. Stensgaard, and F. Besenbacher, 
Appl. Phys. Lett. 77, 708 (2000). 

[3] J. van Ruitenbeek, E. Scheer, and H. Weber, Contacting Individual Molecules Using Me- 
chanically Controllable Break Junctions in: G. Cuniberti, G. Fagas, and K. Richter (Ed.), 
Introducing Molecular Electronics (Springer, Berlin, 2005), pp. 253-274. 

[4] V. Rodrigues, T. Fuhrer, and D. Ugarte, Phys. Rev. Lett. 85, 4124 (2000). 

[5] N. Nilius, T. M. Wallis, and W. Ho, Science 297, 1853 (2002). 

[6] C. L. Kane and M. P. A. Fisher, Phys. Rev. Lett. 68, 1220 (1992). 

[7] I. San and H. Schulz, Phys. Rev. B 59, 3040 (1999). 

[8] G. Cuniberti, M. Sassetti, and B. Kramer, Phys. Rev. B 57, 1515 (1998). 

[9] P. Schmitteckert and F. Evers, Phys. Rev. Lett. 100, 086401 (2008) 
[10] S. Schenk, P. Schwab, M. Dzierzawa, and U. Eckern, Phys. Rev. B 83, 115128 (2011). 
[11] E. Runge and E. K. U. Gross, Phys. Rev. Lett. 52, 997 (1984). 

[12] S. Kurth, G. Stefanucci, C.-O. Almbladh, A. Rubio, and E. K. U. Gross, Phys. Rev. B 72, 
035308 (2005). 

[13] L. Chico and L. M. Falicov, Phys. Rev. B 52, 6640 (1995). 

[14] J. C. Cuevas, A. L. Yeyati, and A. Martm-Rodero, Phys. Rev. Lett. 80, 1066 (1998). 

[15] M. Biittiker, Y. Imry, R. Landauer, and S. Pinhas, Phys. Rev. B 31, 6207 (1985). 

[16] A. R. Rocha, V. M. Garcia-Suarez, S. W. Bailey, C. J. Lambert, J. Ferrer, and S. Sanvito, 

Nature Mater. 4, 335 (2005). 
[17] A. R. Rocha, V. M. Garcia-Suarez, S. W. Bailey, C. J. Lambert, J. Ferrer, and S. Sanvito, 

Phys. Rev. B 73, 085414 (2006). 
[18] J. M. Soler, E. Artacho, J. D. Gale, A. Garcia, J. Junquera, P. Ordejon, and D. Sanchez-Portal, 

J. Phys.: Condens. Matter 14, 2745 (2002). 
[19] Y. Meir and N. S. Wingreen, Phys. Rev. Lett. 68, 2512 (1992). 
[20] A. Chakrabarti, Phys. Lett. A 366, 507 (2007). 

[21] S. H. He, H. U. Baranger, and W. Yang, J. Chem. Phys. 126, 201012 (2007). 

[22] T. Wunderlich, Diploma thesis, Universitat Augsburg (2009), unpublished. 

[23] Note that we use "central region" , "scattering region" , and "scatterer" synomymously. 



12 



